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Abstract 

A  single  collision,  tlma-of'>fiight  mtansion  of  the  VEM  .method  for  the  study  of 
molecule-surface  vibrational  energy  transfer  is  introduced.  The  technique  helps 
election  beMreen  possible  alternative  trial  analytic  forms  of  the  coliisionai  transition 
probability  function.  A  gaussian  form  is  preferred  over  a  boitzmann-exponentiai  form 
for  cyclobutene  isomerization  to  1.3-butadiene  energized  by  collisions  at  a  silica 
surface  at  800  K. 
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MTnooucnoN 

Vibrational  anargy  transfar  during  coliiaions  of  a  thermal  distrioutlon  of  gas 
molecules  with  a  solid  surface  is  a  sublect  of  current  interest.  A  variety  of 

substrate  molecules  have  been  examined  on  various  surfaces  including  fused  silica 

s 

{1.2L  pyrex  glass  C3].  poiycrystalline  iron  (41,  platinum  [5]  and  gold  (61.  Some  of 
the  work  (2.4.51  concerns  measurements  of  vibrational  accomodation  coefficients,  i.e.. 
probes  the  lower  end  of  the  vibrational  energy  distribution  resulting  from  collision. 
The  work  in  this  laboratory  (1.3.61  has  dealt  with  vibrational  excitation  of  substrate 
molecules  to  levels  above  the  unimoiecular  reaction  threshold  Eq  ("*30-65  kcai 
moie~^).  in  this  case.  acts  as  the  absorbing  barrier  postulated  in  the  original 
random  walk  treatment  of  Rubin  and  Shuler  (71. 

The  surfaces  used  in  the  work  cited  above  are  not  'clean*,  being  rapidly 
covered  at  the  substrate  pressures  employed.  In  particular,  they  are  deliberately 

‘seasoned*  at  high  temperatures  in  our  work:  the  seasoned  surface  is  ostensibly 

covered  with  a  graphitic  polymeric  layer  whose  nature  may  be  related  to  those 

described  by  Somorjai.  et  al  (81.  However,  the  nature  of  the  underlying  matrix. 
e.g.  whether  gold  or  silica,  may  continue  to  influence  the  nature  of  the  layer 
formed  by  surface  treatment  and  Its  energy  transfer  properties  (61. 

The  variable  encounter  method  (VEM)  used  in  this  laboratory  has  been 
described  in  detail  (11.  including  the  single  collision  variant  thereof  (31.  in  general, 
cold  'molecules  thermaiized  at  some  lower  temperature  are  permitted  to  collide  with 
a  heated  surface,  called  the  'reactor*,  in  fitting  the  observed  energy  transfer 
results  of  VEM  studies,  various  trial  analytic  forms  of  the  collisionai  transition 
probability  matrix  P  have  been  used  (1.3.91.  Prominent  among  these  have  been 
Exponential  (E).  Qaussian  (Q)  and  Boitzmann-exponentiai  (BE)  forms:  the  elements  of 
the  latter  two  obey  the  following  relations  for  down-transitions  (energy  loss)  by  the 


molecule: 
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p®(AE)  «  C  «xpI-<AE-AE^p)^/2o^]  (1) 

twh«r«  3E  >  Ej  -E|  (Ej  >E| ):  AE^p  Is  the  most-probable  transition:  a  is  the 
standard  deviation;  C  Is  a  normalization  constant; 

p®j  <AE.T^)  «  C  B,  exp(-AE/<AE»  (2) 

where  B|  •  0|  hxp(-E|  /RT^)  is  a  Boltzmann  weighting.  being  the  hot  reactor 
surface  temperature:  <AE>  is  an  average  energy  transfer  parameter:  normalization 
extends  over  all  energies,  and  up-transitlon  probabilities  are  found  by  detailed 
balance,  it  transpires  that  several  transition  probability  models  will  often  fit  a 
particular  range  of  results  although  the  E  form  has  been  found  more  apt  for 
weaker  collisions,  and  the  Q  form  for  stronger  collisions;  the  BE  form  will  obviously 
approach  the  strong  collision  limit  when  <AE>  becomes  very  large. 

In  an  attempt  to  distinguish  better  the  suitability  of  various  models.  Arakawa. 

•  * 

et  al  C91  studied  single  collision  excitation  of  cyclobutene  isomerization  on  a  hot 

seasoned  silica  surface  maintained  at  various  temperatures  and  varied  the  initial 

vibrational  energy,  i.e..  the  initial  (cold)  temperature  of  the  substrate,  from  273 

K  to  600  K.  Although  the  E  model  couid  be  ruled  out.  both  the  G  and  BE  models 

gave  reasonable  fit  to  the  data  taken  at  silica  surface  'reactor*  temperatures  from 

600  K  to  900  K  (see  Fig.  3  of  ref.  9a).  However,  as  may  be  seen  from  Fig.  1 

(for  reactor  temperature  800  K).  these  two  models  make  different  predictions 

concerning  both  the  magnitude  of  the  accomodation  coefficients  a  (a^  <  1: 

BE 

a  <  1).  and  the  form  of  the  population  distribution  above  Eq.  After  a  single 
collision  the  new  population  vector  is  given  by  *  S  wl*hre  is  the  Initial 
Boltzmann  vector  corresponding  to  the  original  (cold)  temperature.  (For  a  strong 
collider,  becomes  the  Boltzmann  distribution  at  the  reactor  (hoO  tempera- 

hire.)  if  P  is  partitioned  at  E^,  i.e..  £  ■  |p  p  |.  then  obviously  the  reaction 

probability  per  collision  P  depends  only  on  P..  i.e..  up-transitlon  probabilities 

C 
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(•rminating  above  Eq.  By  contrast,  maasuramant  of  a  gives  Information  only  about 
At  lower  anargias,  1^^^  follows  the  equilibrium  population  curve  quite 

closely  (hence  «  1),  while  drops  markedly  below  (hence  <  1). 

Above  E^.  declines  more  rapidly  than  Ngg^<  This  difference  in  slope  provides 

s 

a  basis  for  an  experimental  test  of  the  relative  validity  of  the  two  models. 
Unfortunately,  as  may  ba  saen  from  Fig.  2.  this  difference  in  slope  above  E^  Is 

greatly  reduced  at  lower  reactor  temperatures  so  that  in  this  system  It  provides  an 
adequate  basis  for  experiment  only  at  the  highest  temperature  employed  here.  800 
K.  By  changing  the  flight  distance  between  the  hot  surface  and  a  cold  wall,  the 
flight  time  r  during  which  decomposition  may  take  place  can  be  altered.  Since  the 
specific  decomposition  probability  k(E)  varies  with  energy  In  a  range  of  operational 
accessibility  (Fig.  1).  i.e..  k  ^(D  «  0(t).  it  Is  possible  to  probe  the  form  of 

It  may  be  assumed  [3.0.10]  that  accomodation  of  cyclobutene  on  a  cold 
(••2e5K-330K)  surface  is  complete  in  one  collision  so  that  partial  cooling  and 

multiple  reflection  is  not  a  problem. 
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ExpsrIiiMntal 

Cyclobutene  was  obtained  from  Columbia  Chemicals.  >99.9%.  it  was  further 
purified  by  gas  chromatography. 

The  reactor  (Fig.  3)  consisted  of  a  1-i  spherical  pyrex  flask  provided  with  an 

s 

internaiiy-heated  iathe-aiigned  central  silica  cyiinder  finger  having  a  fiat  ground 
fused  silica  end  H  positioned  at  the  center  of  the  fiask.  The  thin  ground  silica  fiat 
(diam.  3  cm)  constituted  the  hot  surface.  The  internai  cylinder  was  described 
previousiy  (3)  and  was  a  thermostatted  iiquid  metal  bath.  B.  The  wall  of  the 
cylinder  was  doubie  Jacketed  and  cooled.  Facing  the  silica  flat  was  an  overlapping 
flat  brass  surface  C  which  was  the  end  of  a  water  cooled  cylinder  that  moved 
smoothly  on  a  screw  S  with  a  no~twist.  no-piay  action  and  pumped  0-ring  seals. 
O.  The  flight  distance  between  fiats  was  changed  from  0.02  cm  to  5.6  cm. 
Temperature  measurements  were  made  with  two  chromei-aiumei  thermocouples  which 
dipped  into  the  liquid  metal  alloy  inside  the  finger.  The  temperature  of  the  hot 
bath  was  varied  from  550-820  K.  The  hot  silica  surface  was  seasoned  with 
cyclobutene  x  10  ^  torr)  at  800  K  for  one  hour.  The  reaction  pressure  was 
0.8-1  X  10  ^  torr.  so  the  system  wes  at  the  second  order  region.  Cyclobutene 
isomerizes  to  1.3-butadiene.  No  by-product  was  found  on  a  FID  chromatographic 
system.  Every  experimental  point  shown  in  Fig.  4  was  the  average  of  three  to  four 
seperate  measurements.  Samples  were  edmitted  to  the  flask  for  a  period  of  time 
that  produced  1-5%  reaction  usually.  Analysis  was  made  on  a  SCOT  squaiane 
column  or  on  a  3/16'  packed  squaiane  column  on  Chromosorb  P. 


RmuHs  and  Calcuiations 

Four  average  flying  distances  were  tested:  0.046.  0.126.  0.32  and  5.6  cm. 
For  tne  former  three  distances,  the  temperature  of  the  cold  brass  surface  was  ••285 
K.  For  the  longest  distance  (the  finger  completely  retracted),  the  wall  of  the  flask 

s 

was  heated  slightly  by  the  hot  surface  and  was  maintained  at  ••320  K.  This  made 

a  slight  change  in  the  number  of  collisions  at  given  pressure  with  the  hot  surface. 

and  In  the  Initial  Boltzmann  distribution.  The  latter  effect  is  slight  and  could 

readily  be  reduced  to  285  K  by  use  of  the  measurements  in  ref.  9a. 

As  before  flO].  the  unimolecular  reaction  obeyed  the  first  order  taw  with 

respect  to  time.  From  the  experimental  rate  constant  k.  the  reaction  probability  per 

collision  P  is  given  by  (Appendix  A).  P.  >  4kV/ca.  where  V  is  the  volume  of  the 
c  c 

reactor:  a  is  the  area  of  the  hot  surface:  c  is  the  average  speed  of  the 
molecules.  A  P^  vs.  T^  plot  is  shown  in  Fig.  4.  The  plot  of  P^  vs.  flight 
distance,  d.  is  given  in  Fig.  5  for  T^  =  800  K. 

P^  values  were  also  calculated  by  computer  simulation. 

Pg  =  J’g^N^E.T)(l-e"‘‘^^^‘^‘‘*'^)dE  /  /QN\E.T)dE 

k(D  is  the  microscopic  rate  constant  according  to  RRKM  theory:  r(d.T)  is  the  flight 

time:  d  is  the  average  flight  distance.  By  fitting  Fig.  4.  the  parameters  <AE>  and 

AE^  p  to  be  used  in  ^  in  order  to  generate  were  deduced  at  several 

temperatures  over  the  range  used.  in  Table  1  are  given  resulting  values  of 

<AE^>.  the  average  down-jump  step,  and  a.  for  the  two  models  under 

consideration  Calculated  plots  of  P^  are  compared  with  experiment  at  800  K  in 

c 

Fig.  5. 

RRKM  values  of  k(E)  shown  in  Fig.  1  were  calculated  with  use  of  activated 


complex  model  A  given  by  Elliot  and  Frey  (11].  Their  complex  C  has  also  been 
used  but  does  not  give  as  good  concordance:  this  is  considered  further  beiow. 
The  frequencies  used  for  the  molecule  and  complex  are  tabulated  in  Appendix  B. 
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DIacuMion 

Reference  to  Table  1  reveals  reasonable  concordance  between  the  values  of 

this  study  and  the  values  of  ref.  9a  for  both  the  average  down-step  size.  <AE^>. 

and  the  vibrational  accomodation  coefficients,  a. 

» 

The  Q  and  BE  models  when  compared  with  experiment  in  the  plot  of  P. 

c 

vs.  T  in  Fig.  5  reveal  better  concordance  for  the  former  model.  Both  calculations 
were  adjusted  by  a  few  percent  so  as  to  fit  the  800  K  data.  The  absolute  values 
of  k(E)  are  sensitive  to  the  choice  of  activated  complex  models,  and  thus  also  is 
the  goodness  of  fit  of  calculated  values  to  experiment.  Use  of  activated  complex  C 
of  Elliot  and  Frey  would  not  give  as  good  absolute  agreement  but.  nevertheless, 
model  Q  still  gives  relatively  better  concordance  than  does  model  BE.  We  conclude 
that  the  gaussian  collislonal  transition  probability  model  Is  the  preferred  one  and 
that,  on  silica,  vibrational  accomodation  coefficient  values  close  to  unity  are  also 
favored. 

The  only  comparable  data  are  those  of  ref.  2  for  octane  on  fused  silica  for 
which  intermediate  values  ("O.S)  were  reported  in  measurements  made  at  lower 
vibrational  energies,  well  below  reaction  threshold  values.  Since  octane  is  a  higher 
boiling  species,  with  a  larger  Lennard-Jones  force  constant  than  cyclobutene,  a 
smaller  a  value  is  not  expected.  The  discrepancy  may  suggest  an  energy  transfer 
model  Intermediate  between  the  two  tested  here.  Our  use  of  models  such  as  E,  Q. 
etc.,  now  and  in  the  past,  has  been  simply  as  pragmatic  analytic  forms  having 
somewhat  different  properties  and  differing  physical  connotations  (121.  There  is  an 
urgent  need  for  the  efforts  of  theoreticians  to  make  the  admittedly  difficult 
connections  with  solid  state  and  molecular  properties  and  potentials  and  with 
eigenstate  densities.  This  need  exists  also  because  of  the  practical  Importance  of 
thermal  high-energy  phenomena. 


Th«  fit  of  calculations  to  experiment  (not  shown)  is  not  as  good  at  700  K 
and  600  K.  although  the  Q  model  continues  to  be  preferred.  As  mentioned  earlier, 
the  operational  distinction  betvreen  the  two  models  decreases,  but  another  difficulty 
is  the  heavier  and  irregular  dependence  of  k(E>  on  the  complex  model  details  at 

s 

lower  temperature  and  energies  where  lc(E)  ceases  to  behave  as  a  smooth 
monotonIc  function  of  E.  This  difficulty  would  be  minimized  in  the  examination  of 
larger  molecules  having  more,  and  more-dispersed,  vibrational  frequencies,  and  with 
a  consequent  wider  range  of  energies  probed  above  E^  at  given  t  ;srature  - 

apart  from  any  fundamental  difficulties  of  the  rate  theory  at  energ  near  the 

reaction  threshold  (131.  Decomposition,  rather  than  isomerization  syste  would  be 

more  favored  In  order  to  keep  values  of  k(E>  up  in  the  experimental.,  accessible 
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Table  I.  Average  energy  transfer  down-step  <AE^>  (cm  S  and  vibrational 

a 

accomodation  coefficient  a  in  the  singie-coiiision  cyclobutene  system. 


ThOO 

Model 

<aeS 

a 

800 

» 

4330  (4100)*^ 

0.95  (0.94)* 

700 

Q 

4940  (4200) 

0.98  (0.97) 

600 

5550  (5640) 

0.96  (0.96) 

800 

4230  (3500) 

0.26  (0.20) 

700 

BE 

5590  (4450) 

0.33  (0.24) 

600 

6580  (6400) 

0.38  (0.33) 

*  For  an  Initial  temperature  of  285  K. 

^  Parenthetic  values  are  average  values  from  ref.  9a  for  initial  temperatures  over 
the  range  273  K  -  600  K. 

^  Values  from  ref.  9a  for  initial  temperature  of  273  K. 


Appendix  A:  The  calculation  of  collision  rate  and  values. 

When  the  flight  distance  is  long,  the  collision  rate  with  the  hot  surface  Is 
obviously  governed  by  the  flask  temperature.  T^.  If  the  flight  distance  Is  only 
several  millimeters  or  less,  the  speed  distribution  of  the  molecules  inside  the 

s 

hot-cold  gap  Is  Influenced  by  the  temperature  of  the  hot  surface.  A  one¬ 
dimensional  approximation  is  suitable  for  our  case.  A  'reduced*  temperature. 

should  ostensibly  be  used  tl41  instead 
of  Tq  to  estimate  the  collision  rate.  However,  the  usual  equation  for  calculating 

the  P  values  from  the  observed  rate  constant  k  is  still  correct  in  this  case, 

c 

When  the  steady  state  is  established.  N^c^  «  '^f  concentration 

of  the  gas  inside  the  gap:  c  Is  the  average  speed  of  the  gas  Inside  the  gap; 

1/2 

Cg>  (8RT^/irM)  :  Nf  is  the  concentration  of  the  gas  outside  the  gap:  Cf  is  the 

1/2 

average  speed  of  the  gas  outside  the  gap.  c^  =  (SRTg/nM) 

Inside  the  gap  one  has  the  usual  equation  for  calculating  P  ; 

P  3  4k^v  /ca:  where  k^  Is  the  'real*  rate  constant  of  this  isomerization  reaction 
c  g  g  g  g 

inside  the  gap:  is  the  volume  of  the  gap:  a  Is  the  area  of  the  hot  surface.  In 

our  experiment  P.  Is  equal  to  or  less  than  10  ^  and  the  composition  of  the  gas 
c 

inside  the  gap  is  the  same  as  the  composition  of  the  gas  outside.  The  quantity 

actually  measured  Is  the  total  rate  constant  k  based  on  the  whole  flask.  Then, 

^  “  VgV^Vf^Vg^'  O'  I'g  -  ®c’ 

41c[V^(T^To)^/^Vg](To/Ti^)^/Vc^a. 

Because  Is  much  bigger  than  V^.  ~  4kV^ /c^  a. 

One  further  addendum.  The  distance  between  hot  and  cold  surfaces,  i.  is 
not  the  average  flight  distance,  d.  Assuming  a  cosine  desorption  iaw.  then 

d  =  r^’^J'^^^lsinedOd^  /J'^’^j’^^^coseslnedOd^  *2i 


n 


Appendix  B:  Vibrational  Assignment  [111 

Molecule  (cm~^)  >  3058.  3126.  9916.  2955(2).  2933.  1566.  1444.  1426.  1288. 
1276.  1210.  1182.  1113.  1100.  1074.  986.  875.  850.  846.  600.  640.  635.  325. 

Activated  Complex  A  (cm~S:  3040.  3030.  3020.  3010.  3000.  2990.  1370. 

s 

1360.  1350.  1340.  1330.  1320.  1310.  1032.  1020.  1010.  1000.  990.  980.  680.  670. 
660.  335. 

Eq  =  32.0  kcai  mol'V  Reaction  path  degeneracy  >  1. 
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Captions 

Fig.  1.  Vibrational  anargy  populations  for  cyciobutana  for  tha  raaction  condition 

(800.  283). - is  tha  tharmaliy  aquiiibratad  Boltzmann  population  at  800  K: 

~  tha  Boltzmann  population  at  tha  tamparatura  285  K:  - .  is  tha 

population  vactor  aftar  a  singla  collision  for  tha  Boltzmann  axponantlai  modal; 

—  •  •  — .  Nq  is  tha  population  vactor  aftar  singla  collision  according  to  tha 
Qaussian  modal.  k(D  is  also  shown  in  this  Figura.  Eq  ■  11.203  cm~\ 

Fig.  2.  Vibrational  anargy  populations  for  cyciobutana  for  tha  raaction  condition 
(800.  285).  Sama  notation  as  in  Fig.  1. 

Fig.  3.  Schamatic  of  apparatus  showing  hot  H  and  cold  C  surfacas:  Intarnal  bath.  B. 
fittad  with  haatar.  stirrar.  tharmocoupia  walls  and  haat  baffla.  coolad  tripia  wall. 
Tha  opposing  watar-coolad  brass  surfaca  C  antarad  via  pumpad  0-ring  saai 
spaca.  O.  on  a  scraw  action.  S.  that  parmittad  7  cm  of  travai.  (3  is  a 
pyrax/quartz  gradad  saai  and  M  a  quartz/stalnlass  staal  transition.  V  ara  spacad 
viawing  ports  for  distanca  calibration. 

Fig.  4.  Tha  raaction  probability  par  collision  P.  vs.  tamparatura. 

c 

□  .A,  o.  ^  rafar.  raspactivaly.  to  tha  avaraga  flight  distanca  5.6.  0.32. 

0.126  and  0.046  cm: - Is  tha  strong  collldar  curva  (E^  «  32.0 

kcal  mola~^). 

Fig.  5.  P  vs.  d  (log  scala)  at  800  K:  0  .  axparlmantai  points:  - .  Qaussian: 

c 

- Boltzmann  axponantlai  calculations. 
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